Does the CMB prefer a leptonic Universe? 
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Abstract. Recent observations of the cosmic microwave background (CMB) at 
smallest angular scales and updated abundances of primordial elements, indicate an 
increase of the energy density and the helium-4 abundance with respect to standard 
c| big bang nucleosynthesis with three neutrino flavour. This calls for a reanalysis of 

the observational bounds on neutrino chemical potentials, which encode the number 
asymmetry between cosmic neutrinos and anti-neutrinos and thus measures the lepton 
asymmetry of the Universe. We compare recent data with a big bang nucleosynthesis 
code, assuming neutrino flavour equilibration via neutrino oscillations before the onset 
of big bang nucleosynthesis. We find a slight preference for negative neutrino chemical 
potentials, which would imply an excess of anti-neutrinos and thus a negative lepton 
number of the Universe. This lepton asymmetry could exceed the baryon asymmetry 
by orders of magnitude. 
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1. Introduction 

Two cornerstones of modern cosmology are the measurements of the abundance of 
primordial light elements and the observation of the cosmic microwave background 
radiation (CMB). Both are described very well by the hot big bang model. The 
abundance of light elements is inferred from the observation of carefully selected 
astrophysical objects, for example extragalactic HII regions to determine the primordial 
helium abundance. The cosmic microwave background is measured for example via 
satellites like the Wilkinson Microwave Anisotropy Probe (WMAP) at large scales and 
with telescopes like the Atacama Cosmology Telescope (ACT) and the South Pole 
Telescope (SPT) at small angular scales. Both types of observation provide comparable 
results for the baryon density of the Universe from very different epochs. 

Recent observations with ACT and SPT allow us for the first time to also estimate 
the cosmic Helium abundance from the CMB. The measurement of light element 
abundances at late times (today) as compared to the epoch of photon decoupling allows 
us to refine the tests of standard cosmology. The recent CMB data from small angular 
scales also allow us to compare an estimate of the number of relativistic degrees of 
freedom at the time of photon decoupling with an estimate of that number at the time 
of big bang nucleosythesis (BBN). 

For the determination of the cosmic abundance of helium-4, the CMB analysis might 
have an advantage over the measurement of extragalactic HII regions, since it is just one 
global dataset and there was no chemical evolution at the time of photon decoupling. 
However, the limits on the primordial abundance of helium are much tighter from stellar 
observations, but in turn the baryon density of the Universe is much better constrained 
by CMB experiments. The best dataset to describe the primordial abundance of light 
elements is a combination of both. 

In this work, we re-investigate the possibility of non-standard big bang 
nucleosynthesis, based on SPT results [1] and the recent reinterpretation of the helium- 
4 and deuterium abundance [2J, [3], [4j [5] . We use stellar observations and CMB data 
to constrain the influence of a possible neutrino or lepton asymmetry. To do so, we 
compare and combine different results for the abundance of primordial light elements 
with theoretical expectations including neutrino chemical potentials. We assume that 
neutrinos are Dirac fermions and that they are relativistic before and at the epoch of 
photon decoupling, i.e. m Ui < 0.1 eV, i = 1,2,3. 

2. Large lepton asymmetry 

Standard big bang nucleosynthesis (SBBN) relies on the baryon to photon density. It 
is commonly defined as the difference of the number density of baryons ni, and anti- 
baryons ni normalized to the number density of the photons n 7 : 775 = {nb — ni)/n^. The 
observed rjt, = O(10~ 10 ) shows a tiny excess of baryons, thus a baryon asymmetry. 

SBBN ignores a possibly large lepton asymmetry, hidden in the three active neutrino 
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flavour. The common model assumption is that sphaleron processes equilibrated the 
total lepton and baryon asymmetry in the very early universe and neutrino oscillations 
result in the equilibration of any lepton flavour asymmetry. Together both assumptions 
result in a tiny, unobservable lepton asymmetry. 

However, the existence of sphaleron processes has not been established by 
experiment so far and numerous theoretical models predict a significantly larger lepton 
(flavour) asymmetry, cf. [SI El IE]- This provides motivation enough to consider a scenario 
with large lepton asymmetry. In previous works we have investigated the effects on the 
cosmic QCD transition and the freeze-out abundance of WIMPS. [9], [10] - The purpose 
of this work is to inspect the consequences for CMB and BBN predictions. 

BBN predictions would then be modified by including neutrino flavour chemical 
potentials fi Vf , with / = e, p, r. For fixed temperature, the introduction of a chemical 
potential increases the energy density of neutrinos. Introducing these additional energies 
leads to a faster expansion of the early universe. We denote the Hubble expansion 
rate with non-vanishing neutrino chemical potentials by H' and H is the Hubble rate 
without lepton asymmetry. The difference is commonly expressed via the expansion 
rate factor S = H'/H = (p'/p) 1 / 2 , with the corresponding energy densities p and p' . 
The difference in the energy densities is the observed extra radiation energy density, 
commonly expressed as additional neutrino flavour in the effective number of neutrinos 




with N v = 3 for the three neutrino flavour / = e, p, r, and corresponding neutrino 
chemical potentials £/ = p Uf /T u at neutrino temperature T v . The expansion rate factor 
becomes S = + (7AA r eff )/43. Assuming neutrino flavour equilibration through 
neutrino oscillations before the start of BBN ensures p e = p^ = p T at T = T B bn [II]- 
Note that the normalization of the additional radiation energy density via neutrino 
chemical potentials does not mean, that the effect might only be due to neutrino 
properties. Any other extension, for example a variation of the gravitational constant, 
might also lead to large neutrino chemical potentials and might also have a compensating 
effects [12]. We will concentrate here only on neutrino asymmetry induced chemical 
potentials. Assuming relativistic neutrinos and a lepton asymmetry much bigger than 
the baryon asymmetry \l\ 3> b, but still \l\ <C 1, one can link the neutrino chemical 
potentials to the lepton asymmetry / [9], 

Zf = jr = \ l Y^ (2) 

where s denotes the entropy density. 

There is also a second effect during BBN, which is due to interactions of electron 
neutrinos with ordinary matter. While all three neutrino flavour chemical potentials 
affect the Hubble rate independently of their sign, the electron neutrino chemical 
potential influences also directly the beta-equilibrium e+p -H- n+u e in the early universe. 
It shifts the proton-to-neutron ratio, depending on the sign of p Ue , and so modifies the 
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primordial abundances of light elements. The two effects can be played against each 
other [13] in a way that one compensates the other one. For further discussions of the 
impact neutrinos have on BBN and CMB c.f. [131 IT5l [TBI fTTj |T8] 

3. Used data and method 

To test the theory of standard big bang nucleosynthesis we have two independent 
possibilities, the analysis of the cosmic microwave background (CMB) and the spectral 
analysis of stellar objects. The baryon to photon density rjio = 10 10 rib/n 7 given by 
the 7 year data analysis of the Wilkinson Microwawave Anisotropy Probe (WMAP 
7yr), 7/10 (WMAP) = 6.23 ± 0.17 [23], is in agreement with the value from the 
observation of primordial deuterium of high redshift, low metallicity quasi stellar objects, 
Vl0 (D) = 6.0 ±0.3 p]. 

The stellar abundance of deuterium is the easiest to trace back to its primordial 
value. It is the lightest bound state and thus burned in all star burning cycles to 
3 He. The observed deuterium abundance at any red shift thus provides a robust 
upper limit on its primordial value. The SBBN prediction from r/ 10 (WMAP) is 
D/H= (2.59 ± 0.15) x 10~ 5 [20]. One seeks to observe young, high redshift and low 
metalicity quasi stellar objects. Nowadays nine objects can be used, and the mean value 
depends on the weighting of the results. [21] find a value of D/H= (2.59 ± 0.15) x 10~ 5 , 
where the value reported in [I] is D/H= (3.05±.22) x 10 -5 , significantly higher than the 
SBBN value. Note that for the latter there is only an overlap with the SBBN prediction 
within their 2 sigma deviations. Neglecting the lower observed abundances would even 
lead to D/H= (3.11 ± 0.21) x 10~ 5 , significantly higher then the SBBN prediction [I]. 

The evolution of the relic abundances of 4 He and 7 Li and the systematic errors in 
their observations are more difficult and introduce possible errors in the determination of 
primordial abundances. The SBBN predictions for 4 He, or equivalently its mass fraction 
Y^(SBBN) = 0.2476 ± 0.0004 [20], are in agreement with observations of low metallicity 
HII regions, Y p = 0.2534 ± 0.0083 [3]. However, the same data set, but using a different 
analyzing method leads to Y p = 0.2565 ± 0.006 [2]. Both values point to higher mean 
values of Y p than predicted by SBBN. 

This is also supported by recent CMB experiments. The analysis of SPT and 
WMAP seven year data with data for the baryonic acoustic oscillation (BAO) and 
measurements of H , with Y p left as a free parameter gives Y p = 0.3 ± 0.03 [1]. 
Considering additionally cluster data, the helium fraction is Y p = 0.288 ± 0.029. Taking 
the SPT and WMAP seven year data alone and leaving both, Y p and N e fi free and 
independent of each other leads to Y p = 0.283 ± 0.045. The tendency of a larger helium 
fraction in CMB data then predicted by SBBN or observed at redshift z ~ 3.5 is also 
supported by WMAP data alone and the Atacama Comology Telescope [22] . 

The direct observation of primordial lithium ( 7 Li/H) differ from the SBBN 
prediction by a factor of 4 to 5. The SBBN prediction is 7 Li/H= (5.071^2) x 10~ 10 
[24] and 7 Li/H= (5.24 ± 0.5) x 10" 10 [20]. All observational data differ significantly 
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from this prediction. [25] found 7 Li/H= (1.58 ± 0.31) x 10~ 10 for a set of halo dwarf 
stars and [26] found 7 Li/H= (1.48 ± 0.41) x 10~ 10 for the abundance in omega centauri. 
This is the so-called Lithium problem, which might have multiple roots, see i.e. [24"1 127] . 
It seems that there are several reasons for a depletion of lithium, however a precise 
understanding of this effect is still missing. We can thus regard the lithium observations 
as robust lower limits on the primordial lithium abundance. 

To constrain further the neutrino chemical potentials, we use the CMB analysis 
of SPT pp. For their analysis they used further cosmological data: the effect of 
baryonic acoustic oscillations (BAO), the WMAP 7 year data analysis (WMAP7), and 
measurements of Hq. To use the constraints from CMB observations on we assume 
for nucleosynthesis calculations that the neutrinos are effectively massless at the time 
of photon decoupling, and that they have no interactions with other particles. Thus 
we can assume £/(Tcmb) = O(^bbn)- We also assume that neutrino flavour is in 
equilibrium before the freeze-out of the neutron-to-proton ration at the begin of the 
BBN epoch. These assumptions lead to one additional parameter for SBBN, a non- 
vanishing {i Vj 

This well motivated extension of the SBBN gives us the possibility to use CMB 
data alone to measure r] w , Y p and iV e ff at the same time and thus to constrain £f 
from the CMB alone. In the following, we apply the extension of the standard model 
with neutrino chemical potentials to predict the relative primordial abundances of light 
elements. We perform our calculations with a full numeric BBN code and compare the 
results to recent reported abundances of light elements, displayed in table [TJ 



Parameter 


Value 


Dataset 


Reference 


100Q b h 2 


2.261±0.054 


SPT+WMAP7-(F p +iV cff ) 


m 




0.2485 ± 0.003 


SBBN (N v = 3) 


m 




0.2534±0.0083 


extragalactic HII regions 


m 




0.283±0.045 


SPT+WMAP7-(F p +iV eff ) 


PQ 


(D/H)xl0 5 


2.59±0.15 


SBBN for tt b h 2 (CMB) 


m 




2.82±0.20 


quasar absorption lines 


m\ 




3.05±0.22 


quasar absorption lines 


m 


( 7 Li/H)xlO~ 10 


k n7+o.7i 

°- v ' -0.62 


SBBN for n b h 2 (CMB) 


m 




1.48±0.41 


Omega Centauri 


[26] 




1.58±0.31 


halo star abundance 


[25 j 



Table 1. Used parameters and values. 



We used a modified version of the PArthENoPE code [29] to calculate abundances 
for varying neutrino chemical potentials £/, equal for all flavour, for a present day 
baryon to photon density 100Q b h 2 = 2.261 ± 0.054 from the combined analysis of SPT 
and WMAP7 with Y p and N c s left free and independent of each other [1]. 
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4. Results 

The results of our calculations are presented in figure 1 and table El Comparing observed 
abundances of 4 He, D/H and 7 Li/H to BBN abundance predictions, we can constrain 
the neutrino chemical potential. 




-1 -0.5 0.5 1 

Figure 1. The red lines represent the calculated primordial abundances of 4 He (Y p ), 
D/H and 7 Li/H as a function of neutrino chemical potential —1 < £/ < 1 and for the 2er 
region of 100n B h 2 = 2.261 ± 0.054 [T]. We assumed £ e = ^ = f T . The blue and black 
dashed lines represent various observational constraints on the observed abundance of 
elements. 

Jf.,1. Helium- 4 

For the 4 He abundance observed in extragalactic HII regions [3] we find an allowed 
2-cr region of —0.09 < < 0.05. The helium abundance inferred form the CMB is 
somewhat higher. The combined analysis of the SPT and WMAP 7 year data pQ in a 
ACDM model with Y p and N e g as free parameter and independent of each other allows 
for -0.46 < & < 0.25. 

Several parameter sets with Y p being a free parameter, but N e g = 3.042 fixed, are 
provided in pQ. Comparing the latter with our results is actually inconsistent, since we 
are interested in N e g > 3.042. Nevertheless, we include these cases for comparison. For 
the analysis of SPT and WMAP 7 year data with N c g = 3.042 fixed but Y p left free, we 
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Element 


Allowed Region 


Dataset 


4 He 


-0.091 < i f < 0.051 


4 He (extra galactic HII) 




-0.461 < £ f < 0.25 


SPT+WMAP7-(y p +iV Gff ) free 




-0.411 < } f < 0.05 


SPT+WMAP7-(Y p ) free 




-0.417 < If < 0.035 


SPT+WMAP7+H0+BAO-(y p ) free 




-0.373 < If < 0.072 


SPT+WMAP 7+HO+B AO +clusters- (Y p ) free 


D/H 


-0.347 < £/ < 0.153 


quasar absorption lines [21] 




-0.524 < if < 0.055 


quasar absorption lines [4] 


7 Li/H 


0.7 < £ f 


Omega Centauri [26] 




0.767 < Cf 


halo star abundance [25] 



Table 2. Constraints on from different observational data. Note that our 
constraints on 7 Li/H arc limited to our calculations < |1|. 

found -0.411 < £ f < 0.05. Adding data of BAO and HO leads to -0.417 < £/ < °- 035 - 
Note, that this constraint on positive neutrino chemical potentials is tighter than that 
from 4 He from extragalactic HII observations. Including also data from galaxy clusters 
in the analysis, one finds —0.373 < £/ < 0.072. Although these tighter constraints 
should not be taken seriously, they demonstrate the potential of upcoming CMB data 
releases to further constrain £f. 

Jf.,2. Deuterium 

For the deuterium abundance we found an interesting difference for the 2-cr overlap of 
the two concurring observational mean values with our calculation (see figure 1). For 
the measurement of D/H of [21], —0.35 < £y < 0.15 is allowed. From the data analysis 
of jl], we find —0.52 < < 0.05. The upper bound is almost identical with the one 
from 4 He abundance of extragalactic HII regions. If we would compare our results to 
the higher D/H= (3.11 ± 0.21) x 10~ 5 reported in g], we find -0.51 < < -0.02. 
This deuterium abundance would exclude positive neutrino chemical potentials within 
2-<t deviation. 

4-3. Lithium 

In the bottom panel of figure 1 we confront our calculation for the 7 Li/H with 
observations. We found an agreement of theory and observed abundance for large 
positive neutrino chemical potentials 0.7 < £/ for J2B] and 0.77 < for the abundance of 
[25] . However, these neutrino chemical potentials that large are excluded by the helium- 
4 and the deuterium data. On the other hand, it is very plausible that lithium has been 
depleted in the course of the galactic chemical evolution [2S]. Thus we conclude, that 
even a large neutrino asymmetry of the Universe, would not solve the lithium problem. 
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4-4- A consistent picture? 

Putting everything together, the most stringent constraint on the neutrino chemical 
potential stems from the helium-4 abundance from extragalactic HII regions, —0.09 < 
£j e < 0.05. This is consistent with all CMB data and with the observed abundance of 
deuterium. The lithium problem remains. However, one should keep in mind, that the 
observations of HII regions may not be as representative for the Universe as the helium-4 
abundance inferred from the CMB. If we rely on the helium-4 abundance from the CMB 
and combining it with the deuterium abundance (taking also systematic uncertainties 
into account), leads to —0.46 < ^ MB ' D < 0.15. Even for this increased interval, the 
lithium problem remains. In both situations the negative values of the are slightly 
preferred. 

5. Conclusion 

Recent CMB data hint towards AN c q > and to a larger abundance of primordial 
helium-4. Here we suggested that both findings could be explained by a lepton 
asymmetry of the Universe, much larger than the baryon asymmetry of the Universe. 
Today this lepton asymmetry would hide in the neutrino background. This scenario 
would have interesting implications for the early Universe, especially at the epochs of 
the cosmic QCD transition [9] and WIMP decoupling [TO] . 

The helium fraction reported by CMB observations results in a negative neutrino 
chemical potential ~ —0.3 and AN e g ~ 0.1. In that case we would live in a Universe 
ruled by anti-neutrinos. From our analysis we conclude that the present abundance of 
light elements and CMB data are not able to rule out £/• = 0, the standard scenario of 
BBN. However, the upcoming releases of new SPT, ACT, WMAP and Planck data, will 
allow us to test the idea of a leptonic Universe. 
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